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SUMMARY 

Thia paper presents a numerical method and the results of a computer program 
for solving an exact, three-dimensional, full-potential equation that models rotat- 
ing and nonrotating inviacid, absolutely irrotational, homentropic flows. Besides 
calculating the flows through an arbitrarily shaped rotor or stator blade row 
mounted on an axisymmetric hub and confined in an axisynsaetric duct, the computer 
program is also capable of analysing flow fields about arbitrarily shaped wing-body 
combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The 
governing equation is solved numerically in a fully conservative form by using an 
artificial time concept, a finite volume technique, rotated type-dependent differ- 
encing, successive line overrelaxation, and sequential boundary-conforming grid 
refinement. An artificial viscosity is added in fully conservative form: and an 
initial guess for the potential field is applied, as determined by a two-dimensional 
cascade analysis. 


INTRODUCTION 

This work is based on the principles used in external transonic aerodynamics 
(Jameson, 1974; Caughey and Jameson, 1977) and represents an extension of the 
author's doctoral research in the field of potential, transonic turbomachinery flows 
( Du 1 i kravich, 1979). 

In axial turbomachinery the pressure ratio across a stage can be substantially 
increased by operating in the transonic speed regime. Analyses of these possibly 
shocked flows should account for their full nonlinearity (Rae, 197b). The simplest 
mathematical model that describes such flows exactly is the so-called full-potential 
equation. This equation can be obtained for rotating or stationary turbomachinery 
geometries from the following analysis. 


ANALYSIS 


The relative coordinate system (x,y,t) is attached for a rotor (fig. 1) which 
rotates at the constant angular speed |il| about the x axis. The free stream ad- 
vances along the same axis at the constant speed U . Let 
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be the relative velocity vector of the fluid with respect to the blade, and let the 
absolute velocity vector be defined as 
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V * v r + (0 x r) (2) 

Then the *um of the inertia, centripetal, Coriolia, and preaaure forcea can be ex- 
pressed in the font ('hi 1952) 

V y x (V x V) - VI - TVS (3) 

Here T it the abaolute atatic temperature, S it the entropy, and rothalpy ia de- 
fined aa 

I - h + (v r • V r - Q 2 r 2 )/2 - H - (3 x r) • V (4) 

where h ia the static enthalpy and H is the total enthalpy. 

The entire flow field can be described with a single potential function 


<P(x,y,z) defined as 

V - VP (5) 

if the condition of irrotationality of the absolute velocity vector is satisfied. 


This requires that VI - 0 and V5 * 0 simultaneously everywhere, or that VI ■ V> . 
Consequently there should be no heat transfer between the fluid and solid surfaces, 
the flow should not separate, and all possible shock waves should be weak. As al- 
ready shown (Caradonna and Isom, 1972; Dulikravich, 1979), the continuity equation 

V • (pV r ) - 0 (6) 

can be written in its full-potential vector operator form 

a 2 ^P - (vp • ?)(VP • VP)/2 + 2(Vp • V)((fl x r) . VP) 


- ((n x r) • v)((n x r) . vp) - o 0 ) 

where a is the local speed of sound. The canonical form (Dulikravich and Caughey, 
1980) of equation (7) is 

aV<P - q 2 <P - 0 (8) 


where s is the relative streamline direction (fig. 1). One can consider this 
second-order, quasi-linear, partial differential equation to be the stecdy-state 
limit of the more general (Carabedian, 1956; Jameson, 1974; Dulikravich, 1979) arti- 
ficial, time-dependent equation expressed in a form suitable for the type-dependent, 
finite difference discretization 
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Here the superscript H designates upstream differencing (used only in the regions 
of locally supersonic relative flow for the purpose of numerically approximating the 
proper domain of dependence), and the superscript E designates central 
differencing; t is the artificial time; and the (n,m) plane is locally orthogonal 
to the relative streamline direction (s coordinate). Equation (9) is iteratively 
solved by using a successive line overrelaxation technique where the iterative 
sweeps through the flow field are considered as successive intervals in the artifi- 
cial time direction. The mixed space-time derivatives are obtained by using a ft- 
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nite difference mixture of old, temporary, and new valuea of <P obtained in the 
iterative sweeping process (Jameson, 1976; Dulikravich, 19 7 9). 

In order of obtain a finite difference evaluation of the derivatives, the flow 
field geometry, the governing equations, and the boundary conditions are transformed 
from the physical space (fig, 1) into a parallelepiped-shaped computational space 
(fig. 2). The periodic flow domain about a single arbitrarily shaped blade mounted 
on an axisymmetric hub and confined in an axisymmetric duct is discretized by a num- 
ber of intermediate axisymmetric surfaces (Dulikravich, 1980a). Every two- 
dimensional periodic surface with the blade intersection contour in the middle is 
then transformed (fig. 3) into a rectangular plane by using conformal mapping 
(Dulikravich, 1979), elliptic polar coordinates, and coordinate stretchings and 
shearings. The uniform grid in the (X,Y) computational plane thus remaps back into 
a periodic body-conforming, quasi-orthogonal grid in the physical space (fig. 4). 

Each distorted, three-dimensional grid cell is mapped into a unit cube by using tri- 
linear, isoparametric local mapping functions (Jameson and Caughey, 1977) of the form 
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b 'iE b p (1 + is P ){1 + V (1 + i V (10) 

P-1 

where the subscript p refers to the value at the cube's corner; that is, 

X « ±1 Y * ±1 Z = ±1 (11) 

P P P 

and b stands for any of the following: x,y,z, or <P. 

In order to use the finite volume technique as defined by Caughey and Jameson 
(1977), equation (8) has to be expressed in terms of the computational (X,Y,Z) coor- 
dinates (fig. 2). If tne geometric transformation is defined as 

[J] * c>(x,y,z)/d(X,Y,Z) and D is a determinant of 1J], the modified contravariant 
components of the relative velocity vector are 


lU r V r M r i T - D[J ] _ 1 Lu r v r w r J T 

If (B ) * l J ] b [J T ] -b , the steady part of equation (9) can be written in the 
following matrix form: 


((a 2 - 9 2 r )/(Dq r ))LU r V r W r ll V>s |T + - 2 (' 
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where 


- a/Dq r )lU r V r W r UV m ,> iB ] T ) - 0 


( 12 ) 


(13) 


^VrV " nz + D[B ^ V xYZ Cr ’ lT 


(14) 


and 


% ‘ < 1/D V lU r V r W r Jt W ! 


(15) 


Equation (b), which is actually equation (8) multiplied by p/a , transform* into 


C ,U r) iX + ( pV r) >Y + ( pW r) >z " 0 (“> 

The boundary conditions in the computational space are the following: On the hub 

surface * 0, on the blade surface ■ 0, and on the duct surface 
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Th« periodic boundary condition applied on the surfaces A'A"B'B" (fig. 1) ia 
4>(X ♦ b,0,Z) « <P(X - b,0,Z), where b <_ s, If the abiolute enthalpy of the fluid 
at upstream infinity ia constant, it followe from equation (4) that the Boat general 
flow kinematic* at upatream infinity can be expreaaed aa + (Cg/r) 

The blade trailing edge ia enforced to be a atagnation line. The finite dia~ 
continuity in the velocity potential at the trailing edge i equal to the circula- 
tion T(r) of the velocity field. This diacontinuity (i.e. , no jua^ in atatic 
preaaure) ia enforced at every point of an arbitrarily ahaped vortex aheet (aurface 
D*E'D"E" in fig. 2) after each iterative aweep through the flow field. 


RESULTS 

A computer program hat been developed on the baaia of a previoua analyaia 
(Dulikravich, 19B0e). Thia program uaea a potential field generated by a two- 
dimenaional caacade analyaia aa an initial gueta for the three-dimenaional potential 
field calculation* (Dulikravich, 1980c). The iterative convergence rate ia accel- 
erated by uaing a three - level consecutive, grid refinement sequence. 

The program was tested for two speed regimes. The first teat case was a sub- 
sonic free rotor of a propeller-type, 100-kW wind turbine, the NASA Mod-0 type 
(Dulikravich, }980b and 1980d). Numerical results showing chordwise distribution of 
the relative Mach number at the tip section are presented in figure 5. Because of 
the lack of available published results for three-dimensional, transonic rotor cal- 
culations, we decided to test an atypical ducted rotor having eight nontvisted, non- 
tapered blades mounted on a doubly infinite cylindrical hub with a hub-tip radius 
ratio r |,/ r t of 0.85. The blades were composed of NACA 0012 airfoil sections 
and had a constant twist (or setting) angle of 4° and a relative chord length 
r t of 0.1. Numerical results obtained at a tip section for |fl| “ 120 rpm, 

(M ) ■ 0. 7458, (p ) * lb N/cn/, and (T ) ■ 325 K are shown in figure 6. 

" ^ - v -Cl) 

Insufficient sharpness of the isentropic shock developing on the suction surface is 
due to the relatively coarse grid used and the low number of iteration cycles. 

CONCLUSIONS 

A general computer program has been developed for fast and acc irate, fully con- 
servative finite volume calculations of the full-potential equation for transonic, 
steady, three-dimensional, potential rotating and nonrotating flows. As a part of 
the program, three-level, refined, boundary-conforming grids are generated for arbi- 
trarily shaped cascades of blades mounted on an ax i syrnmet ric hub and confined in an 
axisymmetric duct. The program uses results of the two-dimensional cascade flow 
calculations as an initial guess for the three-dimensional iterative procedure. 
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Artificial viacoaity vaa added in conaervative font, and thia aaaurea the correct 
strength and position of the captured iaentropic shocks. 


♦National Reaearch Council - National Aeronautica and Space Administration 
Reaearch Associate. 
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